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Abstract 

m 

I In this paper we discuss three methods to calculate energy splitting in cosine potential on a circle, 

Bloch waves, semi-classical approximation and restricted basis approach. While the Bloch wave method 
gives only a qualitative result, with the WKB method we are able to determine its unknown coefficients. 
The numerical approach is most exact and enables us to extract further corrections to previous results. 

s 

1 Introduction 

^- 

In this article we study tunneling in quantum mechanical cosine potential with periodic boundary conditions. 
r*| Hamiltonian of this system is given by the formula 

Oh 

± H=\p 2 + ^-{1-cos{2^ 9 X)) (1) 

with X G (0, g 1 / 2 K) where K is the number of minima and with periodic boundary conditions. For shorter 
notation we put V(x) = (1 — cos(27rx)). Within perturbative calculus vacuum energy of the system 
is degenerate. Each ground state is localized in different minimum. Tunneling is responsible for splitting 
of the energies. In general, this effect cannot be studied analytically. One can give a first approximation 
to the splitting using semiclassical-approximation (or WKB approximation) which was developed by G. 
Wentzel, H. Kramers and L. Brillouin in 1926. In this approach, one finds that the ground energy is shifted 
by a quantity which is a nonperturbative function of coupling constant. This classical solution is called an 
instanton. 

Validity of the instanton calculus is limited to systems with widely separated minima with a large potential 
barrier in between. Nevertheless, it has a vast range of applications to modern field theories. In Yang-Mills 
(*C) theory with SU(2) symmetry group zero energy states are pure gauge (see jU [5]). Such fields can be 

viewed as mappings of SU(2) group into itself and they can be divided into sectors due to topological 
properties. These sectors are labeled by topologically invariant Pontryagin index which takes only integer 
values. Clearly, there is no continuous pure gauge transformation connecting two topological vacua belonging 
to different sectors. There are however configurations with non vanishing field strength which interpolate 
between different topological vacua \n), where n is the Pontryagin index. A special configuration is the 
BPST instanton which is the minimal action path connecting \n) and \n + 1) in Euclidean space. Presence 
of the BPST instanton has dramatic consequences for structure of the real vacuum. The true vacuum of the 
theory is not a single topological vacuum \n) but a superposition of all such vacua \8) — e m9 \n), where 
8 is called vacuum angle. 

Such system can be modeled in one dimensional quantum mechanics by a periodic potential. According 
to Bloch theorem, the energy spectrum consists of continuous bands. Each eigenenergy is labeled by an 
angle 8 and energy states are superpositions of states localized in single minima. Unlike in the quantum 
mechanical case, in Yang-Mills theory only one value of 8 is admissible due to superselection rule. No energy 
bands are present and there is a mass gap between the vacuum and the first excited state. The vacuum angle 
in QCD is responsible for violating CP symmetry. On the other hand, there is no experimental evidence for 
CP breaking which imposes a limit on the angle \8\ < 10 -9 [3]. 

The periodic potential in the weak coupling limit shares many features of the double well potential which 
was extensively studied since 1960s |U [31 |5J 00 El] • It was discovered that there are further corrections 
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to the WKB approximation which can be derived from modified Bohr-Sommerfeld quantization condition 
|10| . They come from from multi-instanton molecules (i.e. a classical path in Euclidean space which is 
composed of instantons that are close to each other) and contribute to the ground energy much weaker than 
a single instanton. Secondly, each instanton molecule contribution (including single instantons) is multiplied 
by a series, which is presumably asymptotic. Moreover, as stated recently by M. Unsal [TTJ , interactions 
between instantons can heal non Borel summability of perturbation series for potentials with degenerate 
global minima. All above mentioned effects take place also case of periodic potential. 

One has to keep in mind that instanton considerations neglect perturbative contributions to energies, 
which are much larger. Secondly, there are higher order corrections to instanton contributions which become 
significant at stronger coupling. It is understandable that there is a need to verify statements concerning 
instantons and see in what regime of the coupling constant the instanton picture is valid. In quantum 
mechanics there is a very efficient approach, called restricted basis method [IS]. It originates from the 
variational Tamm-Dancoff method [13 . In the restricted basis approach one takes basis states \n) with n 
smaller than a certain large cut-off and finds energies in this subspace. A price for taking many states, 
rather than a few trial functions as in original Dancoff paper, is that the calculations have to be performed 
numerically On the other hand, it is very efficient in one dimensional quantum mechanics. Indeed, a good 
convergence of energies with growing cut-off has been observed numerically [14 . Accuracy of this technique 
is limited only by precision of computations and the size of the space. Apart from these limitations, this 
method is exact and provides then a powerful tool for testing WKB approximation. 

Plan of this paper is the following. In Section [2] we present the traditional approach to periodic potential 
with the use of Bloch waves. In Sections [3] we apply the instanton calculus in order to find low energies 
of the system. Results of the semi-classical approximation are consistent with preceding section and give 
more quantitative answer. Section [4] provides the most complete answer. It confirms results obtained by the 
WKB method and show further perturbative corrections. A few coefficients of the series are extracted from 
numerical data. 

2 Bloch waves 

Let us first consider system of infinite size, i.e. K = oo. According to Bloch theorem, the lowest energies of 
the system form a continuous band of width A. They are usually parameterized by an angle: 

E{9) = E- y cos 9, 0e(-7r,7r). (2) 

Corresponding wavefunction is a plane wave modulated by a periodic function ug(x) with period g^ 1 / 2 . 
which is the period of the potential: 

ijjg(x) = exp (iQxg 1 / 2 ^ ug(x). (3) 

Both, E and A depend on specific shape of the potential and are not determined by the Bloch theorem. 

For finite K the wavefunction satisfies periodic boundary condition ipe(0) — ipg^g -1 / 2 ) and only discrete 
values of 9 are allowed. They are: 6 — 2nj/K with j = —K/2 + 1, . . . , K/2 for even K and j = —{K — 
l)/2, . . . , (K — l)/2 for odd K. One can see that energies are doubly degenerate: E(9) = E(—9) for all 9 
except 9 — 0, w. 

Bloch theorem gives a qualitative answer to what are low energies of the system. Still, values of the mean 
energy of the band E and width A need to be found. 

3 WKB approximation 

In this section we will show how lowest energies of the hamiltonian may be obtained in the well known 
instanton calculus. Analogous calculation was given in detail for the double well potential by S. Coleman in 
[15j . A general discussion concerning periodic potentials can be found in |16| . These methods will be applied 
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to the special case of cosine potential with periodic boundary conditions. The transition amplitude from 
minimum x = to x = g~ x / 2 k, k = 0, . . . , K — 1 in Euclidean time T may be expressed as a path integral 

(g- 1/2 k\e- TH \Q) = N ( V[x{T)]e- SE ^\ (4) 



where Af is a normalization factor. The integral is over trajectories which satisfy boundary conditions 
x(-T/2) = 0, x(T/2) = g~ 1/2 k. The Euclidean action S E is 



S E [x(r)} = dr (U(r) 2 + g-'Vig^xir)) ) • (■.) 

J-T/2 \ z 



A trajectory x(t) that minimizes the action for k = 1 is called an instanton. For T — > oo analytical 
solution yields x(t) = \g~ x ^ 2 arctan(e r ) with 5*0 = Se[x(t)] — 2/n 2 g. For k > 1 minimal solutions are 
composed of many such instantons and are called multi-instanton paths. Each instanton of such path 
connects two neighbor minima and begins where its predecessor ended. The integral Q is calculated in 
gaussian approximation around multi-instanton paths and yields 

N J V[x(r)]e- S ^ = -L e - T / 2 £c„^ (J^Tj JC n e'^. (6) 

The coefficient c n ^k is the number of topologically different n-instatnon configurations which satisfy appro- 
priate boundary conditions. Constant K, is 



det 



K 



det'[-£ + V"(g-^x(r))} 



1/2 

= 2. (7) 



Determinants in the above formula are understood as products of all eigenvalues and symbol ' indicates that 
the lowest eigenvalue is omitted. Calculation of JC is postponed to Appendix |A} 

Let us now calculate c n ^. Obviously, c 0j fc = 5o,Jfc and the boundary condition is c n k = c n .o- Let x(t) be 
an n-instanton path ending at fcth minimum. Then the penultimate instanton has to end at minimum k — 1 
or k + 1. Therefore. c n .k = c n -i,k-i + Cn— i,fc+l- Recursive relations can be written as a matrix equation 
c n ,k = Yli RkiCn-i.i with Rki = <5/c-i.i + 8k+ij- Let us notice that eigenvalues of R are = 2 cos(2nj / K) 
with corresponding eigenvector (vj)k = exp(2mjk/K). Then. c n ^ = Ylj a j{ v j)k^j- Coefficients ay are 
determined from initial condition c 0i fc = <5o,fc an d yield oij = 1/K for all j. Finally, 



oo K-l 



(g-^k\e- TH \0) = J—e- T l*Y. E ^"'"i [cos(2nj / K)/cJ ^e^T 

V ^ n=0 j=0 ' V y 

^-e- T / 2 V e 2 ^' fe / A 'exp f Acos{2Kj/K)-^e- 2 l* 2 VT 



(8) 



i=0 

One can use the identity \E) (E\ = 1 to expand the amplitude {g~ x l 2 k\e~ TH |0) as follows: 

(ff-^fcle-^lO) = £ ^ 1/2 k\E) (E|0) e"^. (9) 

By comparison with formula ^ we extract energies and values of wavefunctions at minima 

Ej = \- 4cos(27rj/X)-i=e- 2 ^ 2ff , (10) 



2 y' TT 3 g 
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Note that Ej — EfC-j- Therefore, each energy for j = 1, . . . , L^F^J ^ s degenerate. The lowest energy E$ is 
always non-degenerate and Ek/2 is non-degenerate for even K. 

For infinite K the parameter 9 = 2 ^- becomes continuous. Energies and eigenstates take the same form 
as in Q and ((3): 

Ee = l-4cos(0)-L=e- 2 ^ 2 °, (12) 

(.g- 1 / 2 ^) = Tr^V^. (13) 
To each energy Eg = Ez-k-b correspond two states \6) and \2ir — 0). 



According to [ZJ] there are perturbative corrections to formula (12 1: 



Ej =^a fc g fe -2cos(2^7^)W^ e - S0 ^6 fc5 fc +O(e" 2S Mog( ff )) (14) 

fe=0 fc=0 



with ao = h, bo = 1. 



4 Restricted basis approach 

Yet another technique, restricted basis method can be used to obtain eigenenergies of the Hamiltonian. It is 
performed after |14| . Instead of Fock basis we use plane waves which are more convenient in case of periodic 
potential: 

l/4 

(x\n) = !L =exp(2ning 1/2 x/K). (15) 
v K 



The Hamiltonian is symmetric under TLk group transformation. Let T be the shift operator: T \x) = 
\x + g^ 1 ! 2 )- Its action on the basis vectors is T \n) = exp(— 2irin/K) \n) . Since [H, T] — 0, the Hamiltonian 
H can be diagonalized on each eigensubspace of T separately. There are K eigensubspaces of T, Hj K = 
span{\n)j , n € Z} where \n) - = \ j + nK). Then the Hamiltonian is an infinite tridiagonal matrix 

2ir(j + nK)\ 2 s 1 



j{m\H\n) j = g I — J 5 m>n + g^2^( 2 ^,n - *m,n-l - <*m,n+l)) m,neZ. (16) 

Let us note that the hamiltonian _ff in sector 'Hj.if is the same as in the sector HijjK for any / G N. Therefore, 
the spectrum for K = oo contains all energies from sectors Hj k with p and if being coprime integers. The 
set of lowest energies from all sectors form a dense set in an interval which is the continuous energy band 
for infinite K. 

In order to obtain energies one has to introduce a cutoff \n\ < N and use numerical methods to find 
energies. Because the matrix is sparse, Arnoldi algorithm is very efficient. 

One of the most important issues is rate of convergence. We shall note that in all known cases convergence 
of energies is exponential if the spectrum is discrete and it converges roughly like TV -1 if it is continuous. 
For each finite K spectrum is discrete. For K — oo the spectrum is continuous. However, it is discrete 
in each sector Wj t K separately while there is an infinite number of sectors. Convergence in each sector is 
exponential. Energies as functions of cutoff iV are presented in Fig. ^ 

We are interested in calculating width of the energy band A = Ek/2 ~ Eq, where Ej is the ground energy 
in the sector Hj,K- The smallest value of coupling constant used in our computations was g = 9.13 x 10 -6 . 
Precision of computations is determined by value of A for different g. For g — 9.13 x 1CP 6 it is A = 
2.6 x 10~ 9637 . Needed cutoff was N = 13000. 
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Figure 1: Convergence of energies with growing cutoff N for K = A and g = 10~ 4 . Each line represents four 
energies which are almost degenerate. For larger g convergence is even faster. 



5 Comparison of results 

We will first check the agreement between the Bloch wave approach and WKB method. By comparing ^ 
and ( 12 1 one can read mean energy and width of energy band in 



E=\ (17) 
A = A WKB = 8—=e-^«~3. (INI 




From (13 1 we can see that \0) is an eigenstate of the translation operator T with eigenvalue e 16 . Thus, the 
wavefunction (x\8) has the same form as the Bloch wave djjl. 

Comparison of the WKB and restricted basis methods requires more refined analysis. The semiclassical 
method neglects perturbative corrections to energies, which are seen in the other approach. Therefore, we 
will compare only width of the low energy band A, which is purely nonperturbative. The semiclassical 
approximation is valid for small g so we expect that A num /Ayy KB — 1 — > as g — > 0. This convergence is 
shown in Fig. [2] We can read from (14 1 that the ratio A num / A WKB is a power series: 



A num /A WKB = b k9 k - (19) 

It turns out that numerical results are precise enough to extract several coefficients It appears that 
6fc2 5fc 7r _2fe are integers well within error estimates. Their numerical values are 
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Figure 2: Relative difference between semiclassical and numerical results for splitting of energies tends to 
as g — >• 0. Linear convergence on the log-log plot indicate that there are power corrections to the WKB 
prescription. 
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(20) 



Perturbative corrections can be extracted from any Ej. However, perturbative corrections can be found 
more easily up to higher orders by a Rayleigh-Schrodinger perturbation theory in a selected minimum of the 
potential. 

We also check whether energy dependence on the parameter 9 — 2nj/K agrees with Results are 
presented in Fig. [3] One can see a very nice agreement for g = 0.02. It is violated for larger values of 
coupling constant. 

Structure of wavefunctions is consistent with Bloch waves. Indeed, since the state \Ej) is in the sector 
'Hj l Kt it is an eigenvector of translation operator T corresponding to eigenvalue exp(— 2mj/K). Therefore 
the wavefunction ^pj(x) = (x\Ej) satisfies 

(x + g- 1 / 2 ^) = (z|Tt|£,.) = cxp(2^7^) {x\E 3 ) , (21) 

i.e. it is a periodic function modulated by a plane wave. 
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Figure 3: Dependence of energies on the index j for K = 20, g = 0.02. There are 2 nondegenerate and 
9 doubly degenerate energies. Dots represent numerical data. The continuous line is function E(9) — 
E — ^ cos 9 with 9 = 2'Kj/K, E = E^u, A = E K / 2 — E Q . Agreement is weaker for larger values of g. 



6 Summary 

We have considered three approaches to find low energies for the cosine potential in the small coupling limit 
- Bloch waves, instanton calculus and restricted basis approach. They all gave consistent results, both for 
finite and infinite K. Wavefunctions corresponding to each energy also have the same form. While the Bloch 
theorem gave us only qualitative results, we found values of the mean energy E and band width A with the 
use of instanton calculus. 

Next, we established agreement of the WKB approximation and the restricted basis space method. It 
turned out that values of E and A determined by the former method are only zero order approximation. 
We found perturbative corrections to A with the numerical technique while calculating corrections to E is 
trivial. 

Another interesting issue is summability of aforementioned series. The perturbative series of ground 
energy is known to be asymptotic and thus not summable. Though, its Borel sum can be given a meaning 
when one includes correction due to interactions of instantons. The author addresses this problem in |17j . 



The series b^g k also appears to be asymptotic. From coefficients (20 1 one can estimate the asymptotic 
behavior bf. w —1.1 x 2.8 k kl. Borel sum of this series may be given a meaning after including corrections 
due to interaction of triples of instantons. General scheme of rcsummation for all orders is proposed in 
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A Calculation of the determinant 

We will now show the method of calculating ratio of determinants 



det 



dr 2 



+ 1 



1/2 



det'[-^r + V"(ff-V2i(r))] 



(22) 



I ..-I L, = + W 1 {t) and let ^ } (r) satisfy 



dr 2 



L^f{r)=\4\r), 



^(-T/2) = 0, A^(t) 



= 1. 



(23) 



T=-T/2 



We say that A is an eig envalue of U if ij)f{T/2) = 0. It was shown in [15 that for two bounded functions 
Wi(t) and W 2 (r) it holds that 



(24) 



det 


d 2 

dr 2 ' 


\-Wx(t) 


^{T/2) 


det 


r d 2 

. dr 2 ' 


\-W 2 (t) 


4 2) (T/2) 



Let us take Wi(r) = 1 and W 2 (t) = V" (g~ 1/2 x{T)). One can check that the solution of (|23J for i = 1, A = 
is ^\t) = sinh(r + T/2). The two solutions of the equation L%y{T) — are 



yi(r) = cosh 1 (r), 

y 2 (r) = sinh(r) + tcos1i~ 1 (t). 



(25) 



Function i/j q '(r) is a superposition of those and for large T it reads ipo ( T ) ~ \ eT ^ 2 Vi{ T ) + e T ^ 2 2/ 2 ( T )- 
Green's function for operator L 2 is 

G(r,r') = { -^ X PYM T \ ? rT>r : (26) 
v ' ; \ -W VCr)^') for r<r' v ; 

where W = 2/i(t)?/ 2 (t) — j/i(t)?/ 2 (t) = —2 is the Wronskian. Let Ao be the smallest eigenvalue L 2 . Then, 
-.(2)^,, 



(I: r;;;(r/2) = r;r'(r/2) / T/2 dr'G(T/2,r')A Vi 2 ) (r') « 1 - ^A , 
where we used the fact that tp^ « ^q 2 " 1 to calculate the integral. It follows that Ao ~ 8e~ T . Finally, 

1/2 



(27) 



det 



fC = VA 



■^ + 1 



det [-^ + y"( 3 -i/2 S ( r) )] 



\4 2) (T/2)J 



(28) 
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